Magnetic structure at zigzag edges of graphene bilayer ribbons 
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We study the edge magnetization of bilayer graphene ribbons with zigzag edges. The presence of 
flat edge-state bands at the Fermi energy of undoped bilayer, which gives rise to a strong peak in the 
density of states, makes bilayer ribbons magnetic at the edges even for very small on-site electronic 
repulsion. Working with the Hubbard model in the Hartree Fock approximation we show that the 
magnetic structure in bilayer ribbons with zigzag edges is ferromagnetic along the edge, involving 
sites of the two layers, and antiferromagnetic between opposite edges. It is also shown that this 
magnetic structure is a consequence of the nature of the edge states present in bilayer ribbons with 
zigzag edges. Analogously to the monolayer case, edge site magnetization as large asmK 0.2/j.b 
(per lattice site) even at small on-site Hubbard repulsion U ~ 0.3 eV is realized in nanometer wide 
bilayer ribbons. 

PACS numbers: 73.20.-r, 73.20.At, 73.21.Ac, 73.22,-f, 73.22.Gk, 81.05.Uw 



I. INTRODUCTION 

Graphene, the two dimensional allotrope of carbon, has 
recently been attracting a great deal of attention. Since 
its isolation three years ago [l[ a plethora of unusual 
and interesting properties has been revealed 0, 0, 13]. 
From the point of view of fundamental physics, low- 
energy quasi-particles in graphene behave like massless 
Dirac fermions propagating at an effective velocity of 
light v ~ 10 6 ms -2 . A rather unusual physics is then 
observed, where the half-integer quantum Hall effect is a 
paradigmatic example 0, @|- Graphene is also regarded 
with great expectations from the point of view of tech- 
nological applications. Stability and ballistic transport 
on the submicrometre scale, even at room-temperature, 
make graphene based electronics a promising possibility. 

The possibility of creating stacks of graphene layers 
with the accuracy of a single atomic layer, providing an 
extra dimension to be explored, is another advantage of 
graphene for electronic applications. Of particular inter- 
est to us is the double layer of graphene - the bilayer. 
Bilayer graphene has shown to have unusual electronic 
properties, though unexpectedly dissimilar to those ex- 
hibited by its single layer parent. The new type of integer 
quantum Hall effect observed in bilayer graphene 0,13], 
which is induced by chiral parabolic bands, is an example 
of its uniqueness. From the point of view of applications, 
bilayer graphene is even more promising for some elec- 
tronic devices. It has recently been shown that the band 
structure of bilayer graphene can be controlled externally 
by an applied electric field so that an electronic gap be- 
tween the valence and conduction bands can be tuned 
in a controllable way [1, E3, EH- This makes the bilayer 
graphene the only known semiconductor with a tunable 
energy gap and may open the door for potential applica- 
tions on atomic-scale electronic devices fl2| . 

Among the uncommon features of monolayer graphene 
we find the rather different behavior of the two possi- 
ble (perfect) terminations: zigzag and armchair. While 



zigzag edges support localized states, armchair edges do 
not [13, [l3, [ill- These edge states occur at zero en- 
ergy, the same as the Fermi level of undoped graphene, 
meaning that low energy properties may be substantially 
altered by their presence. The self-doping phenomenon 
[l6j ]. the edge magnetization with consequent gap open- 
ing in graphene nanoribbons [13], and half-metallicity 
[181 ] are examples of edge states driven effects. 

The presence of zero energy edge states at zigzag edges 
of bilayer graphene has recently been confirmed assum- 
ing a first nearest-neighbor tight-binding model [li! ]. Two 
families of edge states has been found to coexist in the 
bilayer: monolayer edge states, with finite amplitude on 
a single plane; and bilayer edge states, with finite ampli- 
tude on both planes, and with an enhanced penetration 
into the bulk. As in single layer graphene, bilayer edge 
states show up in the electronic spectrum as flat bands 
at zero energy - the Fermi energy of undoped bilayer. 
These non-dispersive bands gives rise to a strong peak 
in the density of states right at the Fermi energy, which 
brings about the question of spontaneous magnetic or- 
dering due to electron-electron interactions. 

In the present paper we study the magnetic structure 
of zigzag bilayer graphene ribbons induced by electron- 
electron interactions, which are included through the 
Hubbard model. Working within the Hartree Fock ap- 
proximation we show that due to the presence of edge 
states, which induce a strong peak in the density of states 
at the Fermi energy, zigzag bilayer ribbons show edge 
magnetization even for very small on-site electronic re- 
pulsion. Moreover, it is shown that the spin configura- 
tion is ferromagnetic along the edge, with parallel spins 
occurring on both layers, and antiferromagnetic between 
opposite ribbon edges. Such a magnetic ordering can 
be interpreted as being a consequence of the edge state 
structure in bilayer graphene. 

The paper is organized as follows: in Scc.[ll]wc present 
the model and the mean field decoupling used here; for 
a better interpretation of our results we review briefly in 
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Sec, IIIII the edge states for non-interacting zigzag bilayer 
ribbons; in Sec. |IV] we present and discuss the results of 
this work; we close with conclusions in Sec. E] 

II. MODEL AND MEAN FIELD TREATMENT 

The study of the magnetic structure in AB— stacked 
bilayer graphene given here is based on the ribbon geom- 
etry with zigzag edges shown in Fig. [TJ We use labels 1 
and 2 for the top and the bottom layers, respectively, 
and labels Ai and Bi for each of the two sublattices 
in layer i. Each four-atom unit cell (parallelograms in 
Fig. [JJ has integer indices m (longitudinal) and n (trans- 
verse) such that mai + na.2 is its position vector, where 
ai = a(l,0) and a2 = a(l, — V3)/2 are the basis vec- 
tors and a w 2.46 A is the lattice constant. The simplest 
model one can write to describe non-interacting electrons 
in Ai?-stacked bilayer is the first nearest-neighbor tight- 
binding model given by, 

2 

Htb = ^2 Htba + H±, (1) 

with, 

Htb,% = —t ^ a\ a {m,n)[b l ,a(m,n)+bi,<j(m-l,n) + 

b i>a {m,n- 1)] +h.c, (2) 

H± = —t± a\ t<r (m, n)b2,a(m, n) +h.c, (3) 

m,n,cr 

where aj jCr (m, n) [i>, jCr (m, n)] is the annihilation operator 
for the state in sublattice Ai (Bi), i — 1,2, at position 
(to, n), and spin a = j, J.. The first term on the right hand 
side of Eq. (U) describes in-plane hopping, t ps 2.7 eV, 
while the second term parametrizes the inter-layer cou- 
pling, tj_/t <C 1. In order to examine the magnetic polar- 
ization due to electron-electron interactions we add the 
Hubbard term to Eq. [TJ The total Hamiltonian describ- 
ing the bilayer system reads, 

H = H TB + Hu, (4) 

where Hjj represents the on-site Coulomb interaction, 

2 

i—l m,n 

+ b\ T (m, n)&i. T (to, n)b\ i (to, njb^i (to, n)] , 

(5) 

The Hubbard model is a good starting point to study 
magnetism whenever the density of states at the Fermi 
energy is large enough to produce effective screening of 
the Coulomb interaction. This is true for the clean bi- 
layer, where a finite density of states at the neutrality 
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FIG. 1: Ribbon geometry with zigzag edges for bilayer 
graphene. 

point produces some amount of screening in the sys- 
tem [2fl|. It is certainly the case in the presence of zigzag 
edges, where the density of states peak at the Fermi en- 
ergy implies very effective screening. 

The system Hamiltonian in Eq. §4§ is treated here 
within mean-field theory. In the Hartree Fock approx- 
imation the mean-field version of Eq. ((4]) reads, 

Hmf = Htb + H^ IF , (6) 

with 
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Hy F = U^2 [nAi-a(rn,n)al a (m,n)a iia (rn,n) 

i—l m,n,<7 

+ nBi-a{m, n)bl a (m, n)b itCr (m, n)] , (7) 

where firmer (to, n) is the electronic density for spin a =f 
, I at the site of sublattice r = A, B and layer i = 1, 2 of 
the cell (m, n), The electronic spin densities hTi t a(fn,n) 
have to be determined self-consistently through, 

fiAi.a- (to, n) — ( oj (to, n)a ha (to, n) ) , (8) 

\ ' / MF 

nBi,a(m,n) =(b\ a (m,n)bi^(m,ri)j , (9) 

where the average (• • -) MF is done with the mean-field 
Hamiltonian in Eq. |(6]). Quantum fluctuations, which are 
ignored within mean-field theory, are expected to reduce 
the magnetic moments but not to change significantly the 
overall magnetic structure. As a further approximation 
we assume that the self-consistent solution of Eqs. dH]) 
and dl]) is to independent, i.e., 

nAi,a(m, n) =fiAi,a{n) = — ^ h A i,a(m, n), (10) 

m 

nBt,a{m, n) =n B i^(n) = — ^ n B i,a(m, ri), (11) 

TCI 

where L is the longitudinal ribbon length. We can justify 
this approximation here because we are mainly interested 
on the study of edge magnetization when edge states are 
present, and, as we will see in Sec. IIIII edge states are 
homogeneous along the edge. Note, however, that we 



keep the sublattice index in Eqs. lfTO|) and ifTTj). meaning 
that we can still have in-cell inhomogeneity. 

Without loss of generality we assume that the ribbon 
in Fig. Q] has N unit cells in the transverse cross section 
(y direction) with n € {0, . . . , N—l} } and we use periodic 
boundary conditions along the longitudinal direction (a; 
direction) . Noting the translational invariance of the rib- 
bon along the x direction, and having Eqs. (fT0|l and ifTTj) 
in mind, it is easy to diagonalize Hamiltonian ((6|) with 
respect to the m index just by Fourier transform along 
the longitudinal direction, H = ^2 k Hk, with Hk given 
by, 



H k = H 



TB.k 



MF 
U,k ' 



(12) 



where, 



HrB,k — 
2 

" 1 E E a U fc < + e*)M*> n) + b^(k, 



*JL E a U( fc ' n ) 6 2, ff (fc, n) + h e. , (13) 



and, 



H, 



MF 
U,k 



U ^2 ^2 [h Ai ^ a (n)al a (k, n)a it<T (k, n) 

i—l n,cr 

+ n B i,-tr(n)bl a (k, n)b it „(k, n)] , 



with self-consistent spin densities given by Eqs. {Hi 
and Ijlip . which can be rewritten as, 



(14) 



nAi,a(n) =— (al^(k,n)ai !<r (k, n) 

k 



MF 



MF 



(15) 
(16) 



All conclusions presented in Sec. IIVI regarding the mag- 
netic structure of zigzag bilayer ribbons are drawn by 
solving Eqs. l|12H16p . 



III. EDGE STATES IN THE 
NON-INTERACTING LIMIT 

It is shown in Sec. |IV] that the results for the edge 
magnetization of zigzag bilayer ribbons are a consequence 
of the edge state structure found in this system [19j. In 
this section we briefly review the main features of bilayer 
edge states for U — in Eq. |(4]), i.e., in the absence of 
interactions. 

The band structure of a bilayer ribbon with zigzag 
edges is shown in Fig. [2] (a) for N = 400, obtained by 
numerically solving Eq. lfl3|) . We can see the partly flat 
bands at E = for k in the range 27r/3 < ka < 47r/3, 
corresponding to four edge states, two per edge. The 
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FIG. 2: (Color online) (a) - Energy spectrum for a graphene 
bilayer ribbon with zigzag edges for N = 400. (b) - Zoom 
in of panel (a), (c) - Charge density of the edge states at 
k/2n = 0.36. (d) - The same as in (c) at k/2n = 0.364. The 
interlayer coupling was set to t±/t = 0.2 in all panels. 



zoom shown in Fig. [2] (b) for ka w 27r/3 clearly shows 
that there are four flat bands. 

In order to understand the spatial structure of edge 
states in bilayer graphene we solve the Schrodinger equa- 
tion, H TB ,k \^,k) = E^ k \^,k} 7 for E^ k = 0, where fi 
labels the eigenstate index including spin. First we note 
that Hamiltonian Htbm in Eq. lfl3|) effectively defines a 
ID problem in the transverse direction of the ribbon. It 
is then possible to write any eigenstate k) as a lin- 
ear combination of the site amplitudes along the cross 
section, 



\H,k) = 

2 

EE [ a i( k > n ) \<H,k,n,a) 



■ fti(k,ri) \bi,k,n,cr) 



n i—l 



(17) 



where the four terms per n refer to the four atoms per 
unit cell, to which we associate the one-particle states 
|cj, k,n, a) = c\ a (k, n) |0), with = ai l<T ,b ii<T , spin 
a =T,i) and i = 1,2. To account for the finite width 
of the ribbon we require the following boundary condi- 
tions, 

ai(k,N) = a2(k,N) = pi(k,-l) = p2(k i -l) = 0. (18) 

After solving the Shrodinger equation for zero energy and 
the boundary conditions in Eq. lfT8|) we find four possible 
eigenstates per k, where the only nonzero coefficients for 
each of them are given by [l9| : 



ai(k,n) = 0, a 2 (k,n) = a 2 (k,0)Dke l ~ 



(19) 



ai(k,ri) = ai(k,0)D%e- i ^ n , 
a 2 (k,n) = ~a 1 {k,0)D^ 1 ^e- 



'-(n-l) 



l-D\ 



(20) 
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p l {k,n)=fc{k,N-l)D^e ih ? n ' , 2 (k,n)=O; (21) 
and 

Mk,n) = -/3 2 (k,N -^'-^-fe^'-^fn' - 
p 2 (k,n)=f3 2 (k,N-l)D™'e^ n ' ; 

(22) 

where Du = — 2cos(fca/2) and n = N — n! — 1, with 
n' G {0, . . . , N — 1}. As is easily seen, the coefficients 
in Eqs. i|19H22p give convergent wave functions only if 
2ir/3 < ka < 47r/3, in which case they represent zero 
energy states localized at the surface - edge states - and 
provide an explanation for the four flat zero energy bands 
in Fig. [2] (a) and (b). Equations (fl9)l and (|20)) correspond 
to edge states localized at the top zigzag edge in Fig. [IJ 
and Eqs. (j2Tj) and (|22|) are edge states localized at the 
bottom zigzag edge in the same figure. Note, however, 
that the solutions given by Eqs. (|19H22|) are exact eigen- 
states only for semi infinite systems, where the boundary 
conditions given in Eq. (jTHJ) are fully satisfied. In a fi- 
nite ribbon overlapping of the four edge states leads to 
a slight dispersion and non-degeneracy. Nevertheless, as 
long as the ribbon width is sufficiently large, this effect 
is only important at ka ~ 2n/3 and ka ~ 47r/3 where 
the localization length is large enough for the overlap- 
ping to be appreciable [l5[. For completeness we give 
the normalization constants appearing in Eqs. <|19H22|) . 

\a 2 (k,O)\ 2 = \0i(k,N-l)\ 2 = l-Dl (23) 

Mfc,0)| 2 = \(3 2 (k,N- 1)| 2 = (i (24 ) 

An example of the charge density associated with 
Eq. (|20|) is shown in panels (c) and (d) of Fig. [2] for 
t±/t = 0.2, where the |ai(/c,n)| 2 dependence can also 
be seen as the solution given by Eq. (fl9| for \a 2 (k, n)\ 2 , 
apart from a normalization factor. Of particular inter- 
est to understand the magnetic structure due to interac- 
tion effects is the fact that edge states in zigzag bilayer 
graphene are such that at one edge they live only on 
sublattice A whereas at the opposite edge they live on 
sublattice B. 



IV. RESULTS AND DISCUSSION 

In Fig. [3] the results for the local magnetization m = 
"■Ai.t — ^Ai,i, for i = 1,2, are shown as a function of 
the Hubbard parameter U for different ribbon widths N. 
For each ribbon width we have computed the local mag- 
netization at sites of the A sublattice belonging to cells 
n = 0, n = 1, and right at the middle of the ribbon 
(see Fig. [J). The first conclusion we can draw is that 
sites near the edge get polarized even for very small U, 
while sites in the middle of the ribbon behave like bulk 
bilayer [HI . Another interesting feature shown in Fig. [3] 
is that at the considered edge the magnetization of A2 
sites is larger than that of Al sites, an asymmetry that 
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FIG. 3: Dependence of the magnetization m = riAi,\ — UAi,u 
with i — 1,2, on the interaction parameter U for different 
ribbon widths N. The magnetization was computed at sites 
n = 0, n = 1, and at the middle of the ribbon. Solid lines are 
for the upper layer (i = 1) and dashed lines for the bottom 
layer (i = 2). The result for graphite double sheet (bulk 
bilayer) is also shown. 



vanishes away from the edge. We will come back to this 
below. As regards the B sublattice its magnetization (not 
shown in Fig. [3]) is always similar to the bulk result even 
right at the edge (n = 0). However, when we move to the 
opposite edge, the A and B sublattices change roles: B 
sites at the opposite edge get polarized for very small U 
while A sites show the bulk result. The conclusion then 
is that edge magnetization involving different sublattices 
at opposite edges is showing up in zigzag bilayer ribbons, 
even for very small U. In particular we get m « 0.2/is 
right at the edge for U = O.lt ~ 0.3 eV, similar to what 
is found in graphene [T3 |. 

A better understanding of the edge magnetization is 
achieved by fixing U and plotting the local magnetiza- 
tion to = fiTi^ — nn.i across the ribbon section for 
Ti = Al,A2,Bl,B2. This is done in Fig. g] for differ- 
ent ribbon widths and for a fixed interaction parameter 
U = O.li. As is clearly seen, for such a small interac- 
tion only the edges are polarized. Moreover, the edge 
magnetization is opposite on opposite edges - antiferro- 
magnetic arrangement across the ribbon. Also, we can 
see that at the edge starting with cell n — only sub- 
lattice A has a finite magnetization, whereas at the op- 
posite edge only sublattice B has non-vanishing magne- 
tization. Finally, it is also apparent that at each edge 
the non-zero sublattice magnetization has same sign in 
both layers - ferromagnetic arrangement along the edge. 
These observations are consistent with first-principles 
density-functional calculations of the magnetic structure 
of graphitic fragments (infinite number of layers) [12] and 
bilayer graphene nanoribbons [23l |. 

We have seen in Sec. [ill] that bilayer edge states have 
the following property: at the edge starting with n = 
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FIG. 4: (Color online) Magnetization m — nn^ — nn,i along 
the ribbon cross section for the Vi = Al, A2, Bl, B2. Three 
different ribbon widths were considered: N = 50, 200, 400 
from top to bottom. The interaction parameter was set to 
U = 0.11 



they live only on sublattice A, while at the opposite edge 
they live only on sublattice B, as given by Eqs. (|19H22p . 
The above results for the edge magnetization may there- 
fore be attributed to the polarization of edge states in 
order to reduce on-site Coulomb energy. This interpreta- 
tion also provides an explanation for the layer difference 
in local magnetization. As mentioned before, it can be 
seen in Fig.[3]that the magnetization at A2 sites is higher 
than at Al sites for the edge starting with n — 0. If we 
recall Eqs. (fl9)l and lf20j) for the wave function amplitudes 
at the considered edge we immediately see that while the 
two edge state families contribute to A2 only one has 
finite amplitude at Al sites. The same is true for Bl 
and B2 sites, in agreement with Eqs. (pi) and lj22"|). As 
regards the antiferromagnetic polarization between edge 
states living in opposite edges, it guarantees a ground 
state with zero total magnetization, as it is known to be 
the case for the half-filled Hubbard model [HI- 

Finally we note that edge magnetization gives rise to 
a finite gap at the Fermi level, in complete analogy to 
monolayer graphene [T3|. Half-metallicity has been pre- 
dicted for zigzag single layer ribbons due to the edge mag- 
netization and the presence of a finite gap [l8j]. We ex- 
pect that bilayer ribbons also become half-metallic, with 
an extra switching capability owing to the effect of a per- 
pendicular electric field fulfill]. 
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V. CONCLUSIONS 

We have studied the edge magnetization in bilayer 
graphene ribbons with zigzag edges. The presence of flat 
edge-state bands at the Fermi energy of undoped bilayer, 
which gives rise to a strong peak in the density of states, 
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a consequence of the nature of the edge states present in 
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The experimental observation of edge magnetism in bi- 
layer graphene nanoribbons, and possible application as 
graphene-based magnetic nanostructures, has the ID na- 
ture of the spin polarized state as a major drawback. In 
single layer graphene a crossover temperature T x « 10 K 
has recently been estimated for the magnetic correlations 
at zigzag edges [25] . Below T x the spin correlation length 
grows exponentially with decreasing temperature, while 
above T x it is inversely proportional to the temperature. 
This behavior limits the long-range magnetic order to 
~ 1 nm at 300 K. We have shown, within Hartree Fock, 
that zigzag bilayer nanoribbons have a broken-symmetry 
ground state with a finite spin polarization along the 
edges even at small on-site Hubbard repulsion. The pres- 
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penetration into the bulk, should affect the crossover 
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length. Further work is needed to understand to what 
extent. 
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